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Abstract 

We investigate the effect of an extra fourth quark generation and FCNC mediated Z and Z' 
bosons on the rare decay mode — > (jjii^ . In the standard model, this mode receives only h ^ d 
penguin contributions and therefore, highly suppressed with branching ratio 5 x 10~^. This 
in turn makes this mode a very sensitive probe for new physics. We find that due to the above 
mentioned new physics contributions there is a significant enhancement in its branching ratio. 
Furthermore, the direct CP violation parameter which is identically zero in the SM is found to be 
quite significant. If this mode will be observed in the upcoming LHCb experiment, it will not only 
provide a clear signal of new physics but also can be used to constrain the new physics parameter 
space. 

PACS numbers: 13.25.Hw, 12.60.-1, 12.60.Cn 
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The B decay modes provide valuable insights to critically test the standard model (SM) 
and to look for the possible existence of physics beyond the SM. One way of searching for 
new physics (NP) is by studying the rare decay modes which are induced by flavor changing 
neutral current (FCNC) transitions. In the SM, such rare decays arise at the one-loop level 
and thus the study of the same will provide us an excellent testing ground for NP. 

Over the years, there has been profound interest in the search for physics beyond the 
SM. The observed discrepancy between the measured S^Ks and S^Ks 111 already gave an 
indication of the possible existence of NP in the B (pKs decay amplitude and this has, 
in one way, motivated many to carry out intensive search for NP. Although the presence of 
NP in the 6-sector is not yet firmly established, but there exist several smoking gun signals 
2! which will be verified in the upcoming LHCb experiment or super B factories. Therefore, 
it is interesting to examine as many different rare decay channels as possible to have an 
indication of new physics. 

In this paper, we would like to explore the effect of the extra fourth generation of quarks 
and FCNC mediated Z{Z') boson(s) in the rare decay mode B^ —>■ (j)7i^ , which is a pure 
penguin induced process, mediated by the quark level transition b dss. The interesting 
feature of this process is that it is dominated by the electroweak penguin contributions as 
the QCD penguins are OZI suppressed, and therefore expected to be highly suppressed in 
the SM. It, therefore, serves as a suitable place to search for new physics. At present only 
the upper limit of its branching ratio is known [sl 

By{B- 4>n-) < 0.24 X 10"^ . (1) 



This decay mode has been analyzed both in the SM [4] and in various extensions of it 
where it has been found that in some of these new physics models the branching ratio can 
be enhanced significantly from its corresponding SM value. 

In order to discuss the effect of fourth quark generation and FCNC mediated Z(Z') boson, 

n 

we would first like to present the SM result using the QCD factorization [6|. As the decay 
mode B^ 07r^ proceeds through the quark level transition b — > dss and is a pure penguin 
induced process occurring at the one loop level, the relevant effective Hamiltonian describing 
this process is given by 

„ 10 

-H'Jj = ^K,v;,^^c,{f^n , (2) 
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where p = u,c, Ci{fi)^s are the Wilson coefficients evaluated at the 6-quark mass scale and 
Oj's are the QCD and electroweak penguin operators. 

In QCD factorization 6|, the decay amplitude can be represented in the form 

A{B-{ps) He,Pi)rr-{p2)) = -i'^2m^U{e* ■ pb)F^- J] A,(af - ^af,^^) , (3) 

^ p=u,c 

where \p = VpbV*^, the QCD coefficients Og^g ^^-^ are related to the Wilson coefficients as 



defined in js] and F^^ is the form factor describing B ^ n transition. It should be noted 
that the QCD coefficients contributing to —>■ 07r^ are independent of p = u,c, (i.e., the 
virtual particles in the loop). Therefore, one can also represent the above amplitude using 
CKM unitarity Au + Ac + At = 0, as 

A(B- ^ 07r-) = ^^2m^/^(e* • pB)F_f-(0)At(«3 - ^«3,w) , (4) 

where we have now omited the superscripts on a's. The above amplitude can be simplified 
by replacing 2m^e* ■ Pb ^ The branching ratio thus can be obtained using the formula 

BriB- <Pn-) = 07r-)p , (5) 

IbTcrriB 

where tb is the lifetime of B^ meson. Another possible observable in this decay mode is the 
direct CP violation parameter, defined as 

^ r(i?+ ^ 07r+) - TjB- ^ 07r-) 

Y{B+ (/)7r+) + T{B- (1)71-) ■ ^ ^ 

In order to have non-zero direct CP violation, it is necessary that the corresponding decay 
amplitude should contain at least two interfering contributions with different strong and 
weak phases. Since in the SM, this decay mode does not have two such different contributions 
in its amplitude, the direct CP violation turns out to be identically zero. 

For the numerical evaluation, we use the input parameters as given in the S4 scenario 

n 

of QCD factorization approach [6]. The particle masses and lifetime of the B meson are 
taken from [rf. The value of the form factor at zero recoil is taken as -F_^'^(0)=0.28. The 
value of the CKM matrix elements used are 
\Vcb\ = 42.21 X 10^2, \Vcd\ = 0.2271 and 7 the phase associated with Kb as 70°. With these 
values as input parameters, the branching ratio obtained in the SM is 

Bt^^\B- ^(j)n-)= 4.45 x 10"^ , (7) 



71, iKbl = 3.96 X 10-3, \Vud\ = 0.97383, 



which is quite below the experimental upper limit as given in Eq. (flj). 

Now, in the presence of NP, the transition amplitude (jl]) receives additional contribution 
and can be symbolically represented as 

A^{B~ ^ <Ptx-) = A^^' + = A^^'{1 + r e'^ e-'^^-'*'^) , (8) 

where P is the weak phase of the SM amplitude i.e., we have used Vtd = \Vtd\e~'^^ with 
value f3 = 0.375, is the weak phase associated with the NP amplitude and S is the 
relative strong phase between these two amplitudes. It should be noted that the strong 
phases are generated by the final state interactions (FSI) and at the quark level they arise 
through absorptive parts of the perturbative penguin diagrams. Furthermore, r denotes the 
magnitude of the ratio of NP to SM amplitude. Thus, we obtain the CP averaged branching 
ratio (Br) = [Bt{B~ <pT^^) + Br(i?+ 07r+)]/2, including the new physics contribution 
as 

(Br) =Br^^(l + r2 + 2rcos(/?-0)cos(5) , (9) 

where Br^^ is the SM branching ratio. It can be seen from the above equation that if r is 
sizable, the branching ratio could be significantly enhanced from its SM value in the presence 
of new physics. The direct CP violation parameter ([6]) in the presence of NP becomes 

^ ^ 2r sin 6 sin(0 - (3) . . 

l + r2 + 2rcos5cos(0-/?) ■ ^ ' 



We now consider the effect of a sequential fourth generation of quarks [81]. This model 
is an extension of the SM with the addition of a fourth quark generation. It retains all the 
features of the SM except that it brings into existence the new members denoted by (t', h'). 
The fourth up-type quark {t') like m, c, t quarks contributes in the h ^ d transition at the 
loop level and hence will modify the SM result. The effect of fourth generation of quarks in 
various B decays are extensively studied in the literature 

Due to the additional fourth generation, there will be mixing among the new h' quark 
and the three down type quarks of the SM and the resulting mixing matrix will be a 4 x 4 
matrix. Accordingly the unitarity condition becomes + Ac + Ai + Xf = and thus the 
effective Hamiltonian modifies as 



^eff ^ 



Gp 
V2 



v,,v:,Y,c^o, + yt'iyt',Y.^to. 



(11) 
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TABLE I: Values of the new Wilson coefficients at nih scale where (7"'^"' represents C*' for the fourth 
quark generation model and Cj for the FCNC mediated Z boson model. The phase (p' = {(j) — (3) 
is the relative weak phase between the NP and SM amplitudes. 
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where Cf are the new Wilson coefficients arising due to the t' quark in the loop. The 
values of these Wilson coefficients at the My/ scale can be obtained from the corresponding 
contributions from the t quark by replacing the mass of t quark in the Inami Lim functions 
ll| by t' mass (here we neglect the RG evolution of these coefficients from t' mass scale to 



the weak scale My/). 
group (RG) equation 



hese values can then be evolved to the mb scale using renormalization 



12|, as 



Ci{mb) = U^{mb,Mw,a)C{M, 



w 



(12) 



where C is the 10 x 1 column vector of the Wilson coefficients and f/s is the five flavor 
\0 X 10 evolution matrix. The explicit forms of C{My/) and U^{mb, My/.,a) are given in 
1^ . In Table- 1, we present the values of these new Wilson coefficients at scale for a 
representative set of values for m^/ = 400 GeV. 

After obtaining the values of the new Wilson coefficients at the h quark mass scale, one 
can directly write the decay amplitude analogous to (jlj) , due to the fourth generation quarks 
as 

A'''' = ^2mM^* . pB)F^-mt' (a', - i^^^) . (13) 
where 0^3(3 £;vi/)'s cire the new contributions arising from the t' quark contribution. We 



parameterize the new CKM elements as Xf = rae^'^, where is the new weak phase associated 
with Aj. Furthermore, since the unitarity condition has now become modified the elements 
of the 3x3 upper submatrix of the 4x4 quark mixing matrix will be different from 
the corresponding values of SM CKM matrix elements. Since Vtb and Vtd are not precisely 
known (i.e., not directly extracted from the experimental data, but fitted using the unitarity 
constraint) we will use the lower limits from [7| i.e., \Vtb\ = 0.78 and \Vtd\ = 7.4 x 10~^. 

In order to study the effect of the fourth generation, we need to know the values of the 
new parameters (m^/, r^, 0). Based on an integrated luminosity of 2.3/6^^ CDF collaboration 
[isl l gives the lower bound on mf as mf > 284 GeV. Recently it has been shown that the 
observed pattern of deviations in the CP symmetries of B system can be explained in 
the fourth quark generation model if m^/ > 700 GeV [l^. Therefore, in our analysis we 
consider three representative values for rrit' = 400, 600 and 800 GeV. The value of can be 
obtained from the measured mass difference AM^^ of — B^ system and the corresponding 
expression for ^Mb^ in the presence of fourth quark generation can be found in Ref. 10 1. 
Thus, we obtain the values for different m[, consistent with the unitarity condition of 
4x4 matrix as: 3.8 x 10"^ {mf = 400 GeV), 2.7 x 10"^ {mt> = 600 GeV) and 

~ —2.1 X 10~^ [nit/ = 800 GeV). Using these values, in Figure-1 we show the variation 
of the branching (left panel) and the direct CP asymmetry (right panel) with the new weak 
phase (j) for three different values of rrit'. From the figure one can see that the branching 
ratio is significantly enhanced from its SM value and this enhancement is more pronounced 
for large rrit'. It should also be noted that non-zero direct CP violation in this mode could 
be possible in the presence of an additional generation of quarks. 

Now we consider another extension of the SM, where the fermion sector is enlarged by 
an extra down type singlet quark. Isosinglet quarks appear in many extensions of the SM 
like the low energy limit of the Eq GUT models ISj. The mixing of this singlet type down 
quark with the three SM down type quarks provides a framework to study the deviations 
of the unitarity constraint of the 3x3 CKM matrix. The mixing also induces tree level 
fiavor changing neutral currents, which can thus substantially modify the SM results. In 
this model the Z mediated FCNC interaction is given by 1^ 

>C = :^^^[dLaUap7'dLp]Z, , (14) 
Z COS U\Y 
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FIG. 1: Variation of CP averaged branching ratio (l9|) (in units of 10^®) (left panel) and direct 
CP asymmetry (jlOp (in %) (right panel) with the new weak phase 0, where the solid, dashed and 
dot-dashed lines correspond to m^/ = 400, 600 and 800 GeV respectively. The horizontal line in 
the left panel represents the SM value. 



with 



Yl ^^c.y^P = Sap - VLV,p 



(15) 



i=u,c,t 



where a, (3 are generation indices and U is the neutral current mixing matrix for the down 
quark sector. The non-vanishing component of Uafi will lead to the presence of FCNC 
transitions at the tree level. The implications of the FCNC mediated Z boson effect has 
been extensively studied in the context of h physics [l7, Q, 

Because of the new interactions the effective Hamiltonian describing h dss process is 
given as [18|, 



njZ 

^eff 



■^Vt,V;,[C303 + CjOr + CgO, 



(16) 



where the four-quark operators O3, O7 and Og have the same structure as the SM QCD 
and electroweak penguin operators and the new Wilson coefficients C*j's at the Mz scale are 
given by 

1 Um 



Cr{Mz) 
Cg{Mz) 



6 Vt,v;, 

2 Um 

3 Vt,v;, 

2 Um 



Wi 



3 Vt,v:, 



-(1-sin^M- 



(17) 



7 



These new Wilson coefficients wil' 
renormalization group equation [12 



be evolved from the Mz scale to the rrib scale using 
as described earlier. Because of the RG evolution these 
three Wilson coefficients generate new set of Wilson coefficients Ci{i = 3, ■ ■ ■ , 10) at the 
low energy regime (i.e., at the rrify scale) as presented in Table-1. Thus, one can write the 
new amplitude due to the tree level FCNC mediated Z boson effect in a straight forward 
manner from Eq. (jlj) by replacing a^^s^EW) by a^s^EW), where a's are related to the new 
Wilson coefficients C*j(mb)'s. In order to see the effect of this FCNC mediated Z boson 
effect we have to know the value of the parameter Z — h — d coupling parameter which 
can be explicitly written as Uhd = \Ubd\e^'^ and the allowed range of \Ubd\ is found to be 
(2 X 10"'^ < \Ubd\ < 1-2 X 10"^) [19]. In Figure-2, we present the variation of the CP 
averaged branching ratio with \Ubd\ and (p (left panel) and the direct CP asymmetry 
parameter Acp with (j) (right panel), where we have used sin^ = 0.231. From figure-2 it 
can be seen that the branching ratio could be significantly enhanced and large CP violation 
could be possible in this model. 
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FIG. 2: Variation of the CP averaged branching ratio (in units of 10"^) with \ Ubd\ (in units of 10~^) 
and the new weak phase (p (left panel) and the variation of direct CP asymmetry (llOh (in %) with 
the new weak phase (j) (right panel) where the dashed and soUd lines correspond to \Ubd\ = 10"^ 
and 5 x 10~^. 



Now we consider the effect due to an extra U{1)' gauge boson Z'. The existence of 
extra Z' boson is a feature of many models addressing physics beyond the SM, e.g., models 
based on extended gauge groups characterized by additional U{1) factors [20]. Also the 
new physics models which contain exotic fermions, predict the existence of additional gauge 
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boson. Flavor mixing can be induced at the tree level in the up- type and/or down- type quark 
sector after diagonalizing their mass matrices. Here again as in the Z model, FCNCs due 
to Z' exchange can be induced by mixing among the SM quarks and the exotic quark which 
have different Z' quantum numbers. Here we will consider the model in which the interaction 
between the Z' boson and fermions are flavor nonuniversal for left handed couplings and 
flavor diagonal for right handed couplings. The detailed description of the model can be 



found in Ref. 



21 



22], where it has been shown that such model can successfully explain 



the deviations of S^k and S^ir from S^k and also can explain the B ttK puzzle. The 
search for the extra Z' boson occupies an important place in the experimental programs of 



the Fermilab Tevatron and CERN LHC 



231]. At such hadron colliders heavy neural gauge 



bosons with mass upto around 5 TeV can be produced and detected via two fermion decays 
PVipP) ^ Z' ^ I'^l" (Z = e,fi). 

The effective Hamiltonian describing the transition b dss mediated by the Z' boson is 
given by [2^ 

2 



n 



Z' 

eff 



4Gf 

'7f 



(I 



\g1M2 



^tbVtd 



(18) 



where gi = e/ {sin 9w cos 9w) and B^j denote the left (right) handed effective Z' couplings 
of the quarks i and j at the weak scale. The diagonal elements are real due to the hermiticity 
of the effective Hamiltonian but the off diagonal elements may contain effective weak phase. 
Therefore, both the terms in f|T8l) will have the same weak phase due to Bj^f^. We can 
parameterize these coefficients as 



(I 



\giMz 



(19) 



\giMz'J \VtbV;^j 

where (j)' = (j) — (3, (0 is the weak phase associated with -B^^). 

In order to see the effect of Z' boson, we have to know the values of the and or 
equivalently and B^^^. Assuming only left handed couplings are present, the bound on 
FCNC Z' coupling {B^^,) from B^ — B^ mass difference has been obtained in Ref. 24] as 



y\Re{B^,f\ < 5 x 10-^ y\lm{B!^,y\ < 5 x IQ-^ , 



(20) 



where y = {g' Mz / giMz'Y ■ Generally one expects g' / gi ~ 1, if both the U{1) gauge groups 
have the same origin from some grand unified theories, Mz/Mz' ~ 0.1 for a TeV scale 
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neutral Z' boson, which yields y ~ 10~^. However in Ref. 2j] assuming a small mixing 
between Z — Z' bosons the value of y is taken as ?/ ~ 10^'^. Usiiig y ~ 10^^, one can obtain a 
more stringent bound on < 10^'^. It has been shown in \^ that the mass difference of 
-Bs — Bg mixing can be explained if ~ l^t^^tsl- Similarly, the CP asymmetry anomaly 
in i? ^ (pK^-KK can be resolved if \B^^^B^^^\ ~ l^tfe^^sl- From these two relations one can 
obtain ~ 1. Thus, it is expected that S,l,r ~ 10^^. However, in this analysis we vary 
their values within the range (0.01 — 0.001). 

After having an idea about the magnitudes of these new coefficients which are at the 
Mz scale, we now evolve them to the h scale using renormalization group equation 12|. 
The new Wilson coefficients at the scale are presented in Table-1. Using the values of 
these coefficients at h scale we can analogously obtain the new contribution to the transition 
amplitude as done in the case of Z boson. Now using = = ^, in figure-3, we show 
the variation of the CP averaged branching ratio with ^ and the new weak phase (left 
panel) and the direct CP violation with (right panel). In this case also one can have a 
significant enhancement in the branching ratio for large C,, or in other words for a lighter Z' 
boson. Furthermore, the observation of this mode could in turn help us to constrain the Z' 
mass. 

100 1 1 




FIG. 3: Variation of the CP averaged branching ratio (in units of 10^^) with ^ (in units of 10^'^) 
and the new weak phase (p (left panel) and the variation of direct CP asymmetry (jlOp (in %) with 
the new weak phase (p (right panel) where the dashed and soHd Unes correspond to ^ = 10^^ and 
5 X IQ-^. 



In this paper, we have studied the B ^ 07r decay mode in the standard model and in 



10 



some beyond the standard model scenarios. This is a pure penguin rare decay process and 
proceeds through the quark level transition b — * dss, which occurs at the one loop level and 
is therefore expected to be highly suppressed in the SM. The SM prediction of its branching 
ratio is ~ O{10~^) which is below the experimental upper limit of (9(10^''). We have analysed 
this decay mode in the fourth quark generation model and in the FCNC mediated Z and Z' 
models. In the fourth quark generation model, we find that the branching ratio enhances 
from its SM value, with the increasing m^/ and it can have a value of ~ (9(10"*^). In the 
Z and Z' models, the branching ratio can be significantly enhanced for sizable new physics 
couplings \Ubd\ and ^. In these cases it can reach up to (9(10~^) level but still within the 
experimental upper limit. Furthermore, it is found that large direct CP violation could be 
possible in this decay mode in the presence of above mentioned new physics models. Thus, 
if this mode could be observed in the upcoming LHCb experiment it will provide a clear 
signal of new physics and also can be used to constrain the parameter space of various new 
physics models. However, it would not be possible to distinguish between these new physics 
models considering this mode alone. 
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